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ABSTRACT 

Observations of early afterglow emission from gamma-ray bursts (GRB's) with the Swift satellite 
show steep decay of the X-ray light curve, F^{t) cx with a w 2.5 — 4 at ~ 300 — 500 s after 
the burst trigger. The spectrum in this time interval is consistent with a spectrum Fi, oc with 
/3 « 1. Here, we show that these results can be explained as due to emission from the hot plasma 
"cocoon" associated with the jet, which expands relativistically after the jet has broken through the 
stellar envelope, if a substantial fraction of the cocoon kinetic energy is dissipated at scattering optical 
depths tt ~ 10^ — 10^. This results in the bulk of the cocoon photons being observed at X-ray energies, 
after a delay of few hundreds of seconds relative to the gamma-ray photons from the jet. Multiple 
Compton scattering inside the cocoon causes a spread in the arrival times of the X-ray photons. We 
calculate numerically the observed light curve of photons emerging from the cocoon, and show that it 
exhibits a steep decay, which resembles that observed in many GRB afterglows. During the adiabatic 
expansion that follows the dissipation phase, photons lose energy to the expanding plasma, and as 
a result, the emerging photon energy distribution differs from a black-body spectrum, and can be 
approximated as a power law in the Swift XRT band. Comparison of the numerical results with the 
Swift XRT data of GRB050315 and GRB050421 shows good agreement between the light curves and 
spectra during the initial steep decay phase. 

Subject headings: gamma rays: bursts — gamma rays: theory — plasmas — radiation mechanisms: 
non-thermal — radiative transfer 



1. INTRODUCTION 

The successful launc h of the Swift ^aumi a-my bursts 
(GRB) explorer (e.g., iGehrels et 0(11120041) enabled to 
probe X-ray afterglow emission at early times 10^ s 
to ~ 10^ s) after the burst onset, a time interval which 
was largely unexplored previously. The early afterglow 
observations revealed an unexpectedly steep decline in 
the X-ray light curves of ^ 2/3 of the bursts that did 
not show flaring activity, Fy{t) cx with a k, 2.5 — 4 
following the prompt emission and lasting few tenth - 
hundreds of seconds. The steep decline in the light curve 
is follo wed by a much shallowe r decay, a k. 0.7 a. t later 
times UTaghaferri fi^ nl. 1 120051 IChincaTirii et at I l2005t 
Zhang et al. Il200fit iNousek et al. 11200,51 lO'Brien et al. I 
2006ft . The steep decay segment observed by the 
X-ray telescope (XRT; 'Burrows et al. 2005J usually 
connects to the spectral extrapolation of th e Burst 
Alert Telescope (BAT; Barthel mv et al. Il2005aft p rompt 
emission light curve sni oothly ( Cusumano et al. Il2006t 
IRarthelmv nl ll2005HiHi11 et al \\2n()(t). 

The spectral index obtained by the XRT in the 0.3- 
10 keV range during the steep decay segment of the 
light curve is close t o unity, (x with /9 ~ 1 

ijO'Brien et al. \ l200(t) . This power law index is some- 
what softer than the power law index obtained at higher 
energies by the BAT data, at 15-150 keV . 

A possible explanation to the steep dec line observed 
in many light curves was suggested by iZhang et al I 
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and INousek et al I l)2005D (see also iLiang et al I 

According to their suggestion, the steep de- 
cline is due to the ef fect know n as "the cur vature ef- 
fect" p 'enimor e et al. l [l996: Ku mar fc Panaite scu 200(J 
IDermer ..2004,') . In this effect, radiation arising from high 
angular latitude relative to the viewing direction arrives 
to the observer at late times, following the termination of 
the prompt emission, due to the extra distance it travels. 
In addition, this radiation is observed at flux and peak 
frequency lower than the corresponding in the prompt 
emission radiation, due to the beaming effect. This 
model predicts a connection between the temporal and 
spect ral indices of t he emission (Kuma r & Panaitescu] 
l2000HDCT^ er 2003 IZhang al. li2606|) : in the relation 
/^(t) (X t'^v"^, a and /3 are related via a — 2 + f3. 
However, a detailed a nalysis of 40 XRT light curves by 
lO'Brien et al. I l)2006D found that the spectral and tem- 
poral indices are uncorrelated, in contrast to this model 
prediction. Moreover, when fitting the light curves in a 
model th at allows s tretching and shifting of the trigger 
time, O'B rien et al. (2006) found that 6 out of the 40 
bursts show spectral decay index larger than the allowed 
value of the high latitude emission model prediction. 
In recent years, ther e is increasing evidence t hat long 
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are associated with the deaths of ma ssive stars , 
sumably a rising from core collapse (Woosle^l 
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Levinson fc Eichle r 199 3 ; MacFadycn fc Wopslev l[ 

MacFadve n et'^al 1 1200 it IZhang aD l2003t) fsee also 
Pe'er fc Wiiers 1 120051 for short summary of observa- 



tional evi dence) . In this so called "coUapsar" model 
l|Wooslev I IT99a IPaczvhskilll998(l . the collapse of the 
Fe core of a massive star leads to the formation of a 
central black hole (BH) of several solar masses. A rela- 
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tivistic let is thought to a rise during this short accr etion 
period jMov et oi. ll2?)nnl iMacFadv en et al. ll20ni|) . As 
the hght (compared to the steUar density), relativistic 
jet makes its way out of the progenitor star, its rate of 
advance is slowed down, and most of the energy out- 
put during tha t period is deposited in t o a hot cocoon 
surrounding it l|M eszaro s fc Rees II2001I: iMatzner 
IRamirez-Ruiz et al. 11200^ 

The jet head propagates through the Fe and He 
core of the star, until it emerges from the He core 
into the low den sity H envelope at r^, 10^^ cm 
llMacFadve n & W ooslev IIT99 9; Mes zaros fc Rees 1 120011 
IMatzner Il2003t iWaxman fc Meszaros I I2003D . Following 
the jet emergence, the hot plasma composing the cocoon 
(shocked jet material which accumulates while the jet 
forces its way out, entrained with material swept from 
the star) swiftly escapes from the stellar cavity and accel- 
erates in approximately the same way as the plasma com- 
posing the jet - its energy is co nverted via adiabatic ex- 
pansion into bulk kinetic energy f Meszaro s fc Rees 120011 
[Ramirez-Rui z et al. ■ .200Z') . The entropy per baryon in 

the cocoon, rjc = (Lc/Mc^), is smaller than the entropy 
per baryon in the jet, due to the entrainmcnt of the co- 
coon with material from the star. The cocoon therefore 
accelerates up to an asymptotic Lorentz factor which is 
lower than the Lorentz factor of the jet (but still, rela- 
tivistic), before dissipating its kinetic energy. 

The larger number of baryons in the cocoon compared 
to the jet baryon load implies that if the cocoon ki- 
netic energy is dissipated at similar radii to the dissi- 
pation radius of the jet kinetic ene rgy (e.g., by mag- 
neti c reconnection, as sugges ted by iThompson I l|1994fl 
and lGiannios fc Spruit I ()200,'iD . or alternatively by shock 
waves), than the optical depth for scattering of photons 
inside the cocoon at the dissipation radius is very large, 
Tje ~ 10^ — 10^. As a result, these photons cannot emerge 
from the cocoon immediately, and are diffused out dur- 
ing the adiabatic expansion that follows the dissipation. 
These photons thus suffer a significant time delay com- 
pared to photons that were emitted from the jet as the 
jet kinetic energy was dissipated, which produce the ob- 
served prompt emission. Moreover, the multiple scatter- 
ing of these photons inside the cocoon produces a spread 
in the escape times of these photons, which - as we show 
below - can result in the steep slopes in the early after- 
glow light curves observed by the Swift. 

In this paper, we analyze the emission from the ex- 
panding cocoon, and show that it can provide a natural 
explanation to the observed decline in the early afterglow 
light curves. We calculate in f2|the diffusion time delay of 
the cocoon photons. In Sj^we present numerical results of 
the light curves and spectra. We show that the resulting 
spectra is much flatter than black-body spectrum, due to 
the multiple scattering by the cocoon plasma electrons, 
during the cooling down of the electrons in the adiabatic 
expansion that follows the dissipation. We compare our 
results with the Swift data in 21 before summarizing in 

m 

2. SCALING LAWS FOR THE TIME DELAY AND 
TEMPERATURE OF COCOON PHOTONS 

Following the emergence of the relativistic jet from the 
star (e.g. a WR star), the hot plasma composing the 
cocoon escapes from the jet cavity. As pointed out by 



iLazzati fc Beeelman I ()2005ft . the jet opening angle at 
breakout is < 1°, while the cocoon opening angle is of 
the order of tens of degrees. Following its emergence, the 
cocoon expands in all directions in the dilute stellar wind 
or dilute H-envclope (if any). At this stage, the lateral 
expansion velocity of the cocoon is close to the speed of 
light, and the coc oon expansion becomes isotropic (e.g., 
lAlov eit al ll20?)(il) . We can therefore approximate the 
cocoon geometry in its expansion phase as spherical. 

During the cocoon expansion, which lasts few - few 
tens of seconds (see below), the inner engine that pro- 
duces the jet may still remain active. Therefore, dur- 
ing the cocoon expansion, fresh jet material continuously 
passes through the cocoon, in a similar way to jet prop- 
agation in the core of the star, thus keeping the cocoon 
channel open. During its expansion the cocoon pressure 
decreases, therefore the jet material passing through the 
cocoon is less tightly collimated, resulting in an increase 
of the jet opening ang le as it passes through the c ocoon 
(e.g., see discussion in ILazzati fc Begelman 1 120051 . 

Dissipation of the jet kinetic energy produces the 
prompt 7-ray emission. This dissipation can be caused 
by e.g., internal collisions or magnetic reconnection oc- 
curring at radii larger than that of the leading edge of the 
cocoon (which is less relativistic than the jet, Tco ~ 10 
vs. Tjet ^ 100). Thus, the cocoon does not interfere with 
the usual jet prompt emission production. The presence 
of features in the early X-ray light curve such as X-ray 
bumps or flares, or a shallow X-ray decay phase, should 
also be compatible with the presence of a cocoon. While 
specific mechanisms for flares and shallow decays are not 
considered here, if these are produced, e.g. by late in- 
ternal shocks or continued injection of energy into the 
jet, a requirement is that the funnel drilled by the jet 
inside the cocoon remains open during the prompt and 
early afterglow emission phase. This is plausible e.g. in 
a late internal shock or continued ejection model, where 
material continues to flow out in the jet. This in turn, 
suggests that disturbances and oblique shocks in the ex- 
panding cocoon, caused by the continued jet activity, 
can lead to additional dissipation throughout the cocoon 
volume, in addition to the cocoon interaction with the 
external medium or wind. Photons emitted in the jet, 
whether prompt 7-rays, early flares or shallow decay, are 
expected to escape immediately, provided as usual that 
they are emitted above the jet photosphere. They are 
also beamed into a solid angle smaller than the jet open- 
ing angle. However, photons emitted in the cocoon sur- 
rounding the jet cannot escape directly, due to the high 
optical depth in the cocoon. 

The cocoon photons are therefore observed with a time 
delay compared to the prompt 7-ray photons emitted 
from the jet (although they can arrive before the photons 
produced by late internal shocks or refreshed shocks in 
the jet which might produce flares and shallow decays). 
This time delay relative to the prompt 7-rays has three 
different sources. First, Atg is the "geometrical" time 
delay, representing the time delay of photons originating 
from different regions in the flow compared to photons 
originating from the line of sight toward the observer. 
Second, Atesc is the time delay resulting from the fi- 
nite time it takes the cocoon plasma to emerge from the 
Fe/He core of the star, following the emergence of the jet. 
Third, Atu is the diffusion time delay, which is the time 
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it takes the cocoon photons to escape from the cocoon 
plasma. 

In order to estimate these time delays, we as- 
sume that both the cocoon and the jet emerge 
from the core of the s tar at r^, = lO^^r^ i i era 
llMeszaros fc Rees I 120011 iWaxman fc Meszarosi 120031: 
iMatzner 120031: iLazzati fc Beeelman l2005j) . After the jet 
drills its way through the core of the star, it propagates 
through the H envelope, if any, and escapes, producing 
the "prompt" 7-ray emission at some usual nominal ra- 
dius r.i ~ 10^'^ri_i3 cm, where the jet is already optically 
thin. Relative to some nominal signal emitted from the 
base of the flow, these photons reach the observer sub- 
ject to a geometrical time delay ^^tgjet — ^i/(4r^c) ~ 
Sri^isF^^ ms, due to the finite width and angular size of 
the jet emitting region. 

The cocoon starts its acceleration once it emerges from 
the core of the star, when the jet has broken through. 
We approximate the time it takes the cocoon plasma to 
emerge from the core as Atesc = fi"i,/cs ~ 6r*^ii/o s, 
where Cs = c/\/3 is the speed of sound and / > 1 is an 
unknown numerical factor which accounts for the excess 
of time (above the free expansion time) it takes the co- 
coon to emerge from the core, and is taken here to be 
f — Ifo- To estimate the total energy in the cocoon, 
we start from the Frail et al (2001) value for the typ- 
ical (coUimation corrected) jet energy, Ej ~ 10^"'^ erg, 
with a spread of about one order of magnitude. This 
refers mainly to the electron component, so conserva- 
tively one can multiply by a factor 3 for the proton en- 
ergy in the jet. The cocoon is fed by the total output of 
the jet while the latter expands subrelativistically inside 
the star, and assuming a factor ^ 3 longer jet lifetime in- 
side the star than outside it, the cocoon total energy may 
be taken as Ec ^ W^^Ec.52 erg. The baryon mass in the 
cocoon is also expected to be significantly larger by an 
uncertain value than that in the jet, due to entrainment 
of baryons from the stellar envelope. The entropy per 
baryon in the cocoon plasma is likely to be lower than 
that of the jet, and is parametrized here as rjc = 1077c,i. 
Following its emergence from the core of the star, the 
cocoon expands spherically and accelerates up to a co- 
coon saturation radius, rg^c = i^*Vc = 10^^ r^^nFc,! cm, 
where Tc = rjc = lOFc^i is the asy mptotic value of th e 
cocoon Lorentz factor (e.g.. 'Ramir ez-Ruiz et aLll2002|) . 
Assuming isotropic expansion of the cocoon outside the 
core, the comoving cocoon width is estimated by Arc ~ 
Per* = lOi^r^^iiFci cm. 

The cocoon energy release can be approximated as an 
impulsive release involving a single pulse of width Atesc ■ 
Thus, at a radius r^^c > ^s.c where dissipation of the 
cocoon kinetic energy occurs (e.g. via internal shocks, 
reconnection, etc), the cocoon comoving proton density 
can be estimated as np(c) — Ec/{4:Trrf ^nipC^), where Ec 
is the energy content of the cocoon. For dissipation at a 
radius r^^c below the cocoon photospheric radius where 
the cocoon becomes optically thin, the cocoon optical 
depth to scattering is given by 

T^e,c = Arcnp{c)aT — 350 £'c,52?'*,ii'',~i3rc,i: (1) 

where Ec = 10^^£'c,52 erg. Thus, for dissipation radii 
anywhere between the stellar radius rj^^n and the co- 
coon photosphere at r ~ 10^'* cm, the optical depth to 



scattering inside the cocoon is high. Therefore, photons 
that are produced inside the cocoon suffer, in addition to 
a geometrical time delay as they emerge from the pho- 
tosphere, Atg^c = rpfi/ATcC, also a diffusive time delay 
before they escape. 

The diffusive time delay is estimated as follows. Fol- 
lowing the dissipation phase, the cocoon expands in its 
rest frame at the speed of sound. The electron tempera- 
ture inside the cocoon is mildly relativistic, kT' /nieC^ — 
10~^ (see below), and the cocoon is radiatively domi- 
nated, therefore the cocoon expansion speed is close to 
the asymptotic value v < c/ V3. The rapid expansion im- 
plies that the mean free path of photons inside the cocoon 
increases with (comoving) time as l{t) ~ lo{vt/ Avc)^ , 
where = ArdT^e.c is the mean free path at the end of 
the dissipation phase, and a constant expansion velocity 
V assumed. For a large number of scatterings, r^e,c ^ 1, 
the distance traveled by a photon before it escapes is 
ct = J (dl(t)/dt)dt^ where t is the comoving time from 
photon production to escape. We thus find that the co- 
moving diffusion time is 

iEEE Ai^°-~ (2) 

and the observed diffusive time is shorter by a factor of 
r. A full analytical treatment presented in S^gives a nu- 
merical pre- factor 3^^^ (for v — c^ ). An alternative calcu- 
lation of the time delay can be carried out using the fact 
that the photons emerge from the expanding plasma as 
they reach the photospheric radius, Vph- We show in 
that the two approaches yield similar results. Nonethe- 
less, the advantage of the calculation as presented here 
is that it can be generalized to the calculation of the 
spread in the photon arrival time, which results in the 
steep decay (see 331 below). 

The comoving temperature of the cocoon plasma at 
the dissipation radius is calculated as follows. During 
the cocoon emergence time from the star Atesc, the co- 
coon expands to a radius ^ Aigsc x Cs, thus the cocoon 
occupies a volume Vc = (47r/3)(/?'^)^- As it emerges 
from the core of the star, the cocoon temperature T'{rc) 
is calculated using Ec/Vc — aT'{rcY, where a is Stefan's 
constant. At r < Tg^c, T'{r) oc , and at larger radii. 
T'{r) cx r~^/^. We thus find that at the dissipation ra- 
dius. 

The normalized cocoon temperature at the dissipation 
radius is therefore equal to 

_ ksT' _r pl/4 -1/12 -2/3^-1/3 .-3/4 

m c2 ~ -^c, 52^*41 ^^^,13 ^c,l h 

(4) 

3. LIGHT CURVE AND SPECTRA OF PHOTONS 
EMERGING FROM THE COCOON 

While the details of the dissipation process are un- 
certain, the high optical depth inside the cocoon im- 
plies that photons undergo multiple Compton scatter- 
ing before they escape, hence the resulting spectrum is 
expected to be quasi-thermal irrespective of the details 
of the dissipation process. Following the dissipation, the 
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plasma expands and cools. During the cooling phase, the 
plasma electrons are coupled to the photons via Compton 
scattering. The total energy in the photon component 
cannot exceed the total energy dissipated by the elec- 
trons at the dissipation phase, which is approximated 
by a fraction < 0.33 of the total dissipated energy. 
As shown bv lPe'er fc Waxm an (200^, photons are ex- 
pected to give up to ~ 2/3 of their energy to the plasma 
during the expansion phase that follows the dissipation. 
The plasma energy therefore increases by a factor smaller 
than ^ 20% of its initial energy during the expansion 
phase. 

In addition to Compton scattering, possible sources of 
electrons energy loss during the expansion phase are syn- 
chrotron radiation and bremsstrahlung emission. Syn- 
chrotron losses depend on the uncertain value of the 
magnetic field. An upper limit on this value dur- 
ing the dissipation phase can be found by assuming 
equipartition (a fraction ~ 0.33 of the cocoon in- 
ternal energy density is converted to magnetic field), 
which results in a value B ~ 10^ G. The charac- 
teristic momentum of electrons having Maxwellian dis- 

1/2 

tribution with temperature 9ei <C 1 is 7/3 ~ . 
The cyclo-synchrotron power radiated by these elec- 
trons is therefore Psyn,o = i4,/9)[q^B^{jl3)^]/{mlc^) ~ 
10~^ Bq 9ei,-3 ergs~^, where 9ei = lO~^0ei,-3 is used (see 
cq. 0}. The magnetic field, though, decreases quickly 
during the expansion phase. Denoting by R(t) the char- 
acteristic cocoon radius at time t, the number and energy 
density of electrons decrease as nei{t),Uei{t) cx R{t)~^, 
thus B(t) cx R{t)~^^^ at times t > to , where marks 
the beginning of the expansion phase. For t ^ to, 
R{t) ~ Cs{t — to) where is the speed of sound (see 
below), therefore Psyn{t) ex. t~^. We thus conclude, 
that shortly after the beginning of the expansion, at 
t > {Arc/cs){Psyn,otdyn/mec'^)'^/^ the clcctron energy 
loss time becomes longer than the dynamical time, which 
is comparable to the comoving diffusion time tdyn — 
Atf^-. As a result, synchrotron losses do not affect the 
dynamics of the plasma during the expansion phase. 

A similar conclusion can be drawn for bremsstrahlung 
emission. The rate of energy loss of electron 
by bremsstrahlung radia tion during the diss i pation 
is given by (see, e.g., IR vbicki fc Lightman I Il97fl() 

Pbrem = {27T^ 0ei / 3)^/^2^ qy3hm,c^)neigB ^ 1-3 X 
10"^6'g/!_3nei,i4.5 ergs"\ where n^i = np{c) = 

10^'''^nei^i4.5 cm~'^ is the electrons number density at the 
dissipation, and qb is the averaged Gaunt factor. Since 
the electrons number density decreases quickly with time, 
we conclude that bremsstrahlung losses are insignificant. 
We can therefore approximate the expansion as adia- 
batic. 

During the expansion phase, the expansion velocity of 
a plasma element (in its comoving frame) depends on 
the radial position of the element. The surface of the 
plasma expands at the velocity of sound, while internal 
parts expand at lower velocity, v{r) cx r. The coupling 
of the expansion velocity to the position of a plasma cl- 
ement significantly complicates the problem"*. In order 

^ ISun vacv & Titarchuk (1980) calculated analytically the light 
curve of photons escaping from a steady (non-expanding) plasma 
in the limit r-ye ^ 1, by solving the diffusion equation for pho- 



to calculate the emergence light curve and spectra, we 
therefore use a numerical code, based on Monte-Carlo 
method ^ . 

3.1. The numerical model 

We consider a three-dimensional uniform, homoge- 
neous plasma ball that expands in its rest frame at radial 
velocity v{r) — [r/R{t)\cs, where = c/-\/3 is the speed 
of sound, and R{t) = Rq + Cgt is the ball radius at co- 
moving time t after the dissipation (for convenience, we 
assume that the dissipation phase ends, and the adia- 
batic expansion begins at t = 0. Thus, Rq is the comov- 
ing plasma width at the end of the dissipation phase). 
The plasma is characterized by an initial optical depth 
to scattering T-yeit = 0) = tq. As a result of the expan- 
sion, the number density of particles inside the plasma 
decreases as n^i cx R{t)~^, and the optical depth de- 
creases as Tj^{t) cx R{t)nei{t) cx R{t)~'^ . The particle in- 
side the plasma assume a Maxwellian distribution with 
initial temperature given by equation 0] As discussed 
above, particle cooling during the expansion is insignifi- 
cant compared to adiabatic cooling. Therefore, at t > 0, 
the particles cool adiabatically, hence their temperature 
decreases as 9ei{t) cx R{t)~'^. 

Since during the dissipation phase energy is dissi- 
pated inside the entire comoving width of the cocoon, 
we assume here that photons are injected uniformly in- 
side the ball at < = 0. We consider two types of en- 
ergy distributions of the injected photons: (a) mono- 
energetic injection, with initial (comoving) energy eg = 
{ksT' [r.i) / rrieC^) , where T'{ri) is the electrons tempera- 
ture at the dissipation radius, given by equation and 
(b) power law spectrum of the injected photons above eo, 
with a power law index dn^/ds cx £~^. 

The characteristic energy of synchrotron photons emit- 
ted during the expansion phase is at least seven orders 
of magnitude lower than the plasma energy (e^yn ^ 
10~^ eV for B ~ 10^ G), therefore these photons are ne- 
glected. Similarly, the power radiated in bremsstrahlung 
emission during the expansion is at least three orders 
of magnitude lower than the emitted power during the 
dissipation phase. We therefore neglect in the numer- 
ical calculations additional sources of radiation during 
the expansion phase, and consider Compton (and inverse 
Compton) scattering only during this phase. 

Photons interact with electrons inside the plasma ball 
via Compton scattering (no absorption is considered) . In 
the calculation, the exact differential Compton cross sec- 
tion is used. For each escaping photon, the code traces 
its energy and the delay of the photon escape time com- 
pared to a hypothetical photon that was injected at the 
plasma center at i = and did not suffer any scattering 
before escaping the plasma. 

3.2. Numerical results 

The light curves obtained for three different values 
of the optical depth, r^e = 10^,10^,10'^, and comoving 

tons. This method, however cannot be implemented here because 
of the coupling between the comoving plasma radial velocity and 
the radius of the plasma, v{r) oc r. 

^ We use a versio n of the Monte-Carlo code used in 
IPe'er Waxman 1 120041^ for the calculation of photons energy loss 
during the adiabatic expansion. 
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width Rq — 10^^ cm in all three cases are presented in 
Figure n The cocoon isotropic equivalent dissipated en- 
ergy is Ec.52 — 0.3,1,3 respectively, and we assume a 
redshift z = 1 in a flat universe for calculating the ob- 
served flux. These values correspond to a kinetic energy 
dissipation phase that occurs at radii « 10^'^ cm (see eq. 

The mean time delays of photons in the three pre- 
sented graphs are 8,60,210 s for t^^ = 10\ 10^ 10^ re- 
spectively. These values are in ver y goo d agreement with 
the analytical estimate of eauation lA2l Clearly, the light 
curves at late times decrease according to F^{t) oc t^" 
with a w 3 — 4, irrespective of the exact value of the 
optical depth. Similar ste ep decrease in the light curves 
at late times was found bv lSunvaev fc Titarchuk I l|1980D 
for the case of non-expanding plasma and various ge- 
ometries considered there. This suggests that a steep 
decrease in the light curve is a general property of emis- 
sion from plasma characterized by number of scatter- 
ing higher than few tens (in an expanding plasma, the 
mean number of scattering of a photon before it escapes 
is ri se, o(- Tje, while in a steady plasma n^c. c>c r^^; see, 
e.g. iPe'er fc Waxman ll200|) . 

The emerging photons' spectra are presented in Figure 
13 The photon injection is approximated as monoener- 
getical, with energy equal to the peak energy of the ther- 
mal distribution of the electrons, Eq = Oeirrif-c^ = 500 cV 
(in the plasma frame). The observed energy of the in- 
jected photons is therefore eg^- = eq x Tc/{1 + z) = 
2.5 keV, where Fc = 10 and redshift z — 1 assumed. 
For high value of the optical depth, most of the pho- 
tons lose energy before they escape. As shown in 
IPe'er fc Waxman I pOO^ . for r^e = 100 about 70% of 
the energy lost by the photons is transferred to the bulk 
motion of the expanding plasma, and the remaining is 
transferred to thermal energy of the electrons. This en- 
ergy loss results in the quasi-thermal spectrum obtained 
at low energies. Still, for a uniform injection of photons 
as is considered here, the high energy tail of the emerging 
photon spectra is a power-law like with power law index 
/3 ~ 1 {F^ oc i^^^), in the energy range 0.3—10 keV, which 
is the Swift XRT energy band. The case r^e = 1000 pre- 
sented in figure 12 shows an upper limit on the cocoon 
optical depth that is still consistent with the observa- 
tions. Due to the large number of scattering, the photon 
index in the Swift XRT range is not well defined, and 
the data is marginally consistent with a power law index 
/3 ~ 1 in this range. 

A similar result is obtained when one considers the 
more realistic case of power law photon injection above 
the characteristic energy Eq^', as is expected, e.g., if the 
main radiative mechanism is synchrotron radiation, or 
alternatively, multi ple Compton scatt ering of a photo- 
spheric component l)Pe'er et al. II2005I) . In figure the 
dash-dotted line presents the case of power law injection 
of photons above Eq''', with power law index p = 2. In 
this case as well, the obtained spectral slopes have power 
law index /3 ~ 1 in the Swift XRT energy band. 

4. COMPARISON WITH DATA: THE CASES OF GRB050315 
AND GRB050421 

The steep decline in the XRT li ght curves observed in 
many bursts ifO'B rien et al. ''20061 is similar to the steep 
decline in the light curves obtained by the Monte Carlo 
simulation. We show in Figures lHl^ two fits of the data of 
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Fig. 1. — Light curves of photons emerging from an expanding co- 
coon, characterized by initial optical depths Tje,c = 10, 100, 1000. 
Initial comoving radius Rq = 10^-^ cm, comoving temperature 
9^1 = 10""^ injected photon energy eo = 10~^ (in unites of 
nieC^), characteristic Lorentz factor Fc = 10 and redshift 2 = 1 
are assumed. The dissipated cocoon isotropic equivalent energy 
Ec,52 = 0.3,1,3 respectively is considered. Photons are injected 
uniformly inside the sphere. Similar light curves are obtained 
for photon injection with power law energy distribution above eo- 
Each of the presented curves is produced by 10^ Monto-Carlo runs. 
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Fig. 2. — Resulting spectra of the emerging photons for the cases 
of T^e.c = 100, 1000. Solid lines present injection of mono energetic 
photons at energy eo = 10~^ (in units of rrieC^), and dash-dotted 
line presents injection of photons having power law energy distri- 
bution with power law index p = 2 above eo , as might be expected 
from synchrotron emission above the peak and cooling frequencies. 
Thick lines are for T^e = 100, and the thin line presents the case 
of Tje = 1000. All other parameters are the same as in figure 
For T-.,e = 100, the spectral slope in the XRT range (0.3 - lOkeV) 
is close to /3 = 1, while for r-ye = 1000 a single spectral slope is 
marginally consistent with this result. 



GRB050315 and GRB050421 obtained from P. Obrien by 
the numerical model results. These two bursts did not 
show significant flaring activity, and had a good time 
coverage during the steep decay phase of the emission. 

In fitting the data we assumed a cocoon optical depth 
at the dissipation radius T.ye,c = 100. This value is con- 
sistent with the assumption that the dissipation radius 
is > 10" cm. The red shift of GRB050315 is z = 1.949 
l|Vaughan et al. II2006D . and the redshift of GRB050421 
is unknown, thus is taken here to be z = 1. The val- 
ues of the cocoon Lorentz factors during the cocoon dis- 
sipation phase Fc were determined using the observed 
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time delays at the beginning of the decay. The resulting 
Lorentz factors, Tc — 30,8 in the two cases are of the 
same order of magnitude, Tc ~ 10. The value of the dis- 
sipated cocoon energy required to fit the data depends 
on the assumption about the injected photon distribu- 
tion. In the case of GRB050315, an isotropic equivalent 
cocoon energy £'c,52 = 1-5 (£'c,52 = 2.5) is required in 
the monoenergetic (power law injection) scenario, while 
in the case of GRB050421 cocoon energy £^c,52 = 0.03 
{Ec,b2 = 0.1) is required in these two scenarios. The 
values found are consistent with the order-of-magnitude 
estimates presented in f|21based on theoretical arguments. 

The data of GRB050421 (figure El shows clear distinc- 
tion between the early shallow decay observed at the first 
30 s, which is attributed to prompt emission, and 
the rise of the flux followed by a steep fall at ~ 100 s 
which is attributed in our model to emission from the 
cocoon. This time separation is less significant in the 
case of GRB050315 presented in figure |Sl where a flare 
at ~ 25 s precedes the cocoon emission, which peaks at 
~ 30 s. We further discuss this time separation issue in 
SjS] below. 

Fitting of the XRT spectra in these two cases are 
presented in Figures [3 El In the cases presented, we 
use the same parameters of the fits presented in fig- 
ures I3I4I with a re- normalization of the flux to the co- 
coon emission time. We assume initial normalized elec- 
trons temperature O^i — and injected photon energy 
Eo = OgimeC^ — 500 eV in the monoenergetic injection 
case. In the power law photon injection assumption, a 
power law index p — 2 above £o in the photon energy 
distribution is assumed. 

In fitting the data of GRB050315, the numerical re- 
sults were shifted by a factor of 1.2 (3) in the mono- 
energetic (power law) injection cases. A physical origin 
of these shifts is that the temperature of the cocoon is 
different than the assumed numerical value of 6ei = 10~^ 
by equivalent factors. We can thus summarize that the 
comoving temperature of the cocoon is del ~ 8.3 x 10^^ 
(3.3 X 10~^) for the monoenergetic (power law) injections 
assumptions used. Similarly, the comoving temperature 
of the plasma for the case of GRB050421 is O^i = lO'^ 
(5 X 10~^) for monoenergetic (power law) photon injec- 
tion cases. 

For GRB050315, the spectra of the BAT data at late 
time s is fitted with t he XRT data with a single power 
law l|Vaughan et al . 2006). This is consistent with the 
numerical results of the light curve presented in Figure 
|S1 for the more realistic scenario of power law injection 
of photons. BAT data of GRB050421 exists only at the 
first ^ 100 s, before the steep decay phase of the light 
curve. We thus conclude that the BAT photons originate 
from the jet, and not from the cocoon plasma. 

The two fits presented here are only representative 
cases. Similar fits can be obtained for many light curves 
in O'Brien et al. (2006) sample, which show power law 
indices during the decay phase in the range a ~ 2 — 4. 
The similarities found between the numerical model light 
curves for different values of the optical depth imply that 
the data can be fitted with various values of the free pa- 
rameters Tye.c and Tc- However, the demand that Tc is 
of the order of few tens, and not larger than ~ 100, and 
the demand that the spectrum differs from a black-body 
restricts the optical depth in both cases considered here 
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FlG. 4. — Fitting of the decline part in the XRT data of 
GRB050421 with the results of the numerical simulation light 
curve. Here, F;, = 8, -Ec,52 = 0.1 and z = 1 were considered, 
and all other numerical free parameters have the same values as in 
figure 13 



to be r^e,c ~ 100. 

The results presented here thus show that using values 
of the optical depth and cocoon temperature close to the 
"canonical" values considered in fj21(see eqs. HEl): emis- 
sion from the cocoon can account for the observed steep 
decay spectra. The injected photons are assumed in all 
cases to have monoenergetic distribution at energy simi- 
lar to the peak energy of the thermal electrons inside the 
cocoon, or to have a power law distribution above this 
energy. These photons therefore originate by dissipation 
mechanism inside the cocoon, without any assumptions 
on photons produced by the prompt emission. For as- 
sumed cocoon optical depth, by fitting the model to the 
data one can determine the values of the cocoon Lorentz 
factor Tc, the electrons temperature 9ci and the dissi- 
pated cocoon energy, Ec- 

5. SUMMARY AND DISCUSSION 

In this paper, we have considered emission from an ex- 
panding cocoon as a possible explanation to the steep de- 
cay in the Swift XRT light curve observed in many bursts 
after few tens - few hundreds of seconds. We showed in |21 
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Fig. 5.— Fitting of the XRT spectrum of GRB050315 with 
the results of the numerical simulation light curves. Parameters 
used in the fitting are r-ye^c = 100, Ro = lO^'^ cm, Fc = 30 and 
z = 1.949. The cocoon energy and comoving temperature are 
Ec,52 = 1.5, Osi = 8.3 X 10~* for the mono-energetic (solid line), 
^c,52 = 2.5, 0^1 = 3.3 X 10~* for the power law (dash dotted line) 
photon injection assumptions. 
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Fig. 6.— Fitting of the XRT spectrum of GRB050421 with the 
results of the numerical simulation light curves. Parameters used 
in the fitting are T-ye,c = 100, Rq = 10^^ cm, Fc = 8 and z = 1. 
The cocoon energy and comoving temperature are -Ec,52 = 0.03, 
9^1 = 10~^ for the mono-energetic (solid line), Ec,52 = 0.1, S^; = 
5 X 10"^ for the power law (dash dotted line) photon injection 
assumptions. 



that due to the high optical depth in the cocoon, photons 
produced inside the cocoon at similar radii to the esti- 
mated dissipation radius of the kinetic energy of the jet, 
are observed after a diffusion time delay of few hundreds 
of seconds. We then presented in f|31our numerical code, 
and showed the emergent photons light curve and spec- 
tra. We showed that the emergent photons light curves 
and spectra are similar to the decaying part of the light 
curves and the XRT spectra observed in GRB050315 and 
GRB050421, and argued that due to the similarities in 
many of the SwiftXET light curves, the numerical results 
fit the data in many bursts. 

Due to the complexity of the problem, we did not ob- 
tain a simple analytical formulae for the light curves. 
However, we showed that the numerical results can be 
approximated as a power law over a limited time inter- 
val and energy range, which is similar to the energy range 
observed by the Swift XRT. Due to the large number of 



scattering the photons undergo before they escape, the 
resulting spectrum in the XRT band is by large inde- 
pendent on the initial photons spectrum, which by itself 
depends on the details of the dissipation process. We 
therefore conclude that observations of the spectrum dur- 
ing the decay phase are not expected to unveil much of 
the mystery of the dissipation mechanism, if the cocoon 
model is correct. 

We found that the shape of the light curves are similar 
for different values of the optical depth (figure QJ. Light 
curves similar to the light curves presented h ere were 
also obtained bv ISunvaev fc Titarchuk I ljl98CH) . for the 
case of non-expanding plasma. These facts suggest that 
the shape of the light curve may be common to many 
astrophysical phenomena in which photons undergo tens 
to hundreds of scattering before they escape the plasma 
in which they are produced. This may also explain the 
similarities found in the shape of the light curves in ma ny 
different bursts in the sample of lO'Brien et al. 1 lj2006tl . 

These similarities also imply that the data could be fit- 
ted with various values of the free parameters - r^e c, Tc 
and the comoving temperature O^i- We showed in f^two 
fits for the light curves and energy spectra of two differ- 
ent bursts, assuming that photons are produced at the 
dissipation inside the cocoon, without additional prompt 
emission photon source. We found that the fitted val- 
ues of the cocoon Lorentz factor r^, the cocoon energy 
Ec and the electrons temperature O^i are close to their 
"canonical" values calculated in Sj21 based on theoretical 
arguments. Restrictions on the values of these param- 
eters are obtained from the physical demands that Fc 
is of the order of few tens, and not larger than ~ 100, 
and that the plasma temperature should not exceed the 
theoretical values predicted in ^ 

The sample of Swift-XRT light curves presented in 
Figure 3 of O'Brien et al. (2006) shows significant va- 
riety In many cases (e.g., GRB050319, GRB050412, 
GRB050422, GRB050713A, GRB050803 and many 
more), there is a clear distinction between the early 
(prompt) emission phase, and the steep decay phase. 
Following an early shallow decline, the flux rises before 
entering the steep decay phase. This pattern can be at- 
tributed to emission from the cocoon (see figure^. In 
other cases (e.g., GRB050713B, GRB050814 and GRB 
050819) a steep fall in the flux follows the end of the 
prompt emission, without a significant rise preceding the 
decay. In other cases, a shallow decay of the flux is ob- 
served over hundreds - thousands of seconds. 

The model presented here suggests that (at least part 
of) these observations are the result of emission from the 
cocoon. The extension of the prompt emission to tens 
- hundreds of seconds implies that in many cases the 
emission from the cocoon is partially obscured by the 
prompt emission. This may explain the lack of rise in the 
flux in many cases which show steep decay. Obviously, 
the model presented here cannot account for the entire 
variety of the light curves observed. Nonetheless, the 
facts that rise in the flux was observed prior to the steep 
decay in several cases, and the steepness of the decay - 
a « 2.5 — 4 in many cases, are difhcult to account for in 
models that do not consider the existence of the cocoon. 

In addition to the many cases of long bursts which show 
a steep decline in the XRT light curve, a similar behavior 
was observed in one case of a short burst - GRB050724. 
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While the nature of short bursts progenitors is stiU not 
fully understood, it is most likely that these bursts do 
not originate from the explosion of massive stars, hence 
the cocoon model cannot be similarly motivated in this 
case. The decay observed in GRB050724 is so far unique, 
and was not observed in other short bursts. The present 
explanation, which applies to long bursts, does not pre- 
clude a high-latitude emission effect interpretation in 
some bursts, e.g., in short bursts, or in some fraction 
of long bursts. 

Our model assumes an unspecified dissipation mecha- 
nism that converts cocoon kinetic energy into radiation 
at radii < 10^^ cm, below the radius where the cocoon 
becomes optically thin to scattering. This radius could 
be as low as the stellar surface r ^ 10^^ cm, but could be 
anywhere up to the photosphere. Such a dissipation may 
be ca used by magnetic reconnection ( Giannios & Spruitl 
12005ft. interaction with the expanding envelope of the 
star ijMeszaros R,ees ll20?)]l) or extended central engine 
activity and refreshed shocks in the jet that could in- 
teract laterally with the cocoon region, and may derive 
oblique shocks into the cocoon. As presented in sim- 
ilar time delays would in principle occur if the cocoon 
photons were advected from the core of the star, with- 
out any dissipation of the cocoon kinetic energy at larger 
radii. However, in this latter case the very high optical 



depth at ~ 10^^ cm would imply that the emergent spec- 
trum is very close to black-body, which for the bursts 
discussed here is not consistent with the observations. 

In contrast to the internal collision model, which pre- 
dicts multiple dissipation phases resulting in long dura- 
tion prompt emission, the cocoon model predicts only 
a single dissipation phase of the cocoon plasma. This 
model therefore disfavors a second phase of steep decay. 
The flares that are observed in many bursts thus are 
predicted to have a different source, e.g., extended cen- 
tral engine activity. Nonetheless, as discussed above, if 
the dissipation of the jet kinetic energy occurs at small 
enough radii where the optical depth is higher than a few 
tens, similar light curves are expected. 
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APPENDIX 

THE DIFFUSIVE TIME DELAY IN AN EXPANDING PLASMA 

Assuming that following the dissipation phase the cocoon plasma expands at velocity close to its asymptotic velocity, 
Cg = c/\/3 the mean free path of photons increases with comoving time as l{t) = lo{l + vt/ Ar)^, where Iq = Arc/r^e^c 
is the mean free path at the end of the dissipation phase, and constant expansion velocity v = Cg assumed. For a large 
number of scattering, T^e,c 2> 1, the distance traveled by a photon before it escapes is 



cAtT = 



At; 



dl , , 
—dt = la 
dt 
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where A<™' is the comoving photon diffusion time. We thus find that in the observed frame, the diffusion time is given 

by 



Atd = 



At'" 



^ V Fa 

ocn 17I/2 3/2 -3/2t^1/2 
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The diffusive time delay calculated above is of the order of few hundreds of seconds for the assumed values of the 
free parameters. This time delay is much longer than the escape time delay of the cocoon plasma from the core of 
the star Atesc, and it is longer than the geometrical time delay Atgjet of the prompt radiation from the jet, but it is 
comparable to the geometrical time delay from the photosphere of the cocoon, Atg^c, as long as the optical depth to 
scattering inside the cocoon is high. 



ALTERNATIVE APPROACHES FOR THE CALCULATION OF THE DIFFUSION TIME DELAY 

Photons that are produced in an expanding plasma in regions of high optical depth (i.e., at dissipation radii r^) <C r^^) 
emerge from it as the plasma reaches the photospheric radius, rp^. The observed time delay of these photons compared 
to photons that did not suffer any scattering can be calculated in two alternative ways: (a) calculation of the diffusion 
time delay in an expanding plasma, as was done here; and (b) since these photons are advected outward with the flow 
from the emission radius ro to the photospheric radius, rph, the time delay should be equal to the "geometrical" time 
delay of photons emerging from rph- Here, we show that the two different approaches lead to the same result, provided 
the optical depth at the dissipation radius is high and that the plasma expands at velocity close to c. 

In order to include a variety of possible dissipation mechanisms (such as, e.g., magnetic reconnection), we assume 
that the comoving width of the plasma at the dissipation radius is arbitrary, and denoted by Arjj. At larger radii 
r > rz) the dissipated region is advected with the flow, therefore the width of this region increases with radius, Ar cx r. 
For a single (not continuous) energy release, the number density of protons at radius r is n,p{r) = E / (AiTr^mpC^), thus 



the optical depth at radius r is 



EutAtd 
ATTrnr'^rnpC^ 



(Bl) 



Radiation form cocoon 
where Ar(r) = ArD(r/r]j). The photospheric radius is given by rph — r{T — 1), oi 



1/2 



^^'^ \47rrDmpc2 ' ' ^^'^'^ 



The diffusive time delay is approximated by 



A.Ob. _ Vr(r = rj)Arr, r^^ Ar^ Arjrph) , . 

We thus find that the diffus ive time delay is similar to the time delay resulting from the finite width of the plasma at 
the photospheric radius fsee lWaxman 11199^ . In addition, photons that are emitted at the photospheric radius from a 
point not on the line of sight suffer a time delay ~ Vph/T^c, which is thus the minimum observed time delay. For the 
specific case of dissipation by internal shock waves, the comoving width is Ar = r/T, the time delay due to the finite 
width of the plasma is similar to the time delay due to emission from positions off the line of sight, and we obtain 
the familiar result Af^' ~ Vph/V^c. However, the above analysis shows that in the more general case of arbitrary 
plasma width Ar > r/F, the diffusive time delay is equal to the time delay due to the finite width of the plasma, 
which determines the observed time delay. 
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